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FIG. 1: Energy-band diagram of the structure. 

Both the hole and electron reservoirs are assumed to 
be in thermal equilibrium. For the physical phenomena 
we are interested in, the fermi level of the p(n)-side hole 
(electron) is slightly lower (higher) than the hole (elec- 



2 



tron) subband in the dot. After a hole is injected into 
the hole subband in the quantum dot, the n-side electron 
can tunnel into the exciton level because of the Coulomb 
interaction between the electron and hole. Thus, we may 
assume three dot states 



(n D ) t - (n D ) = -ig I dt'{(p) t , - (pt) } - 2T D f dt'(n D ) t , 
Jo Jo 



10} = \0,h) 
\U) = \e,h) 

\D) = |0,0) (1) 

, where |0, h) means there is one hole in the quantum dot, 
e, h) is the exciton state, and |0, 0} represents the ground 
state with no hole and electron in the quantum dot. One 
might argue that one can not neglect the state |e, 0) for 
real device since the tunable variable is the applied volt- 
age. This can be resolved by fabricating a thicker barrier 
on the electron side so that there is little chance for an 
electron to tunnel in advance. We can now define the 
dot-operators % = \U) (U\ , n D = \D) (D\,p= \U) (D\ , 

su = |0) (U\ , sh = |0) (D\. The total hamiltonian H of 
the system consists of three parts: the dot hamiltonian, 
the photon bath, and the electron (hole) reservoirs: 

H = Hq + Ht + Hy 

A A 

-Ho = £u n U + + tip + H res 

k 

Hp = y^^ujkb\bk 

k 

H v = ^(VqC^Su + W q 4»D + cc.) 
q 

q q 

In above equation, bk is the photon operator, gDk is the 
dipole coupling strength, X = ^2 k Dkb\ , and c q and d q 
denote the electron operators in the left ad right reser- 
voirs, respectively. Here, g is a constant with a unit of 
the tunneling rate. The couplings to the electron and 
hole reservoirs are given by the standard tunnel hamil- 
tonian Hy, where Vq and Wq couple the channels q of 
the electron and the hole reservoirs. If the couplings to 
the electron and the hole reservoirs are weak, then it is 
reasonable to assume that the standard Born-Markov ap- 
proximation with respect to these couplings is valid. In 
this case, one can derive a master equation from the exact 
time-evolution of the system. The equations of motion 
can be expressed as (cp. [22]) 



(P)t-(p)t = -T D j* tfe**-V (X t xlptf)) v 
-ig f* dife ie ^{(n v X t xl) v 
-(n A D Xl,X t ) t ) 

(h-ik = -T D j\t>e-^{${t>)X v Xi) v 
+ig j\t'e- i ^- t '\(n u X t ,xt) v 
-(n A D xtx t ,) t ). (3) 

, where IV - 2ir £ q V^S( £u - e%) , T D = 
2lT J2q W q S ( £ D - £q), and e = Eu - e D is the energy 
gap of the quantum dot exciton. Here, p(t') = pe let X t i, 
and X t i denotes the time evolution of X with H p . The 

A 

expectation value (p^) t describes the decay of an initial 
polarization of the system and plays no role for the sta- 
tionary current. Therefore, we shall assume the initial 

expectation value of p vanishes at time t = 0. 

As can be seen from Eq. (3), there are terms like 

nxjX t X\^i which contain products of dot operators 

and photon operators. If we arc interested in small cou- 
pling parameters here, a decoupling of the reduced den- 
sity matrix p(t') can be written as 

p(t') « p° ph Tr ph p(t')- (4) 
By using the above equation, we obtain 

Tr(p(t')nuX t xl) « (4} t , (*t4) ( 5 ) 

and correspondingly the other products of operators can 
be obtained also. For spontaneous emission, the photon 
bath is assumed to be in equilibrium. The expectation 

value (XtX\^ = C(t—t') is a function of the time inter- 
val only. We can now define the Laplace transformation 
for real z, 



(nu) t ~ (nu) = -ig / dt'{(p) t , - (p f ) t ,} 



+2Tu [ dt' 
Jo 



(1 - (nu) t , - (n D ) t ,) 



C e (z) 



nu(z) 



Jo 

-r 

Jo 



dte~ zt e iet C{t) 



die zt (nu) t etc., z > 



(6) 
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and transform the whole equations of motion into z- 
space, 

nu(z) 

n D (z) 

p(z) 
p*(z) 

These equations can then be solved algebraically. The 
tunnel current / can be defined as the change of the oc- 
cupation of ri(7 and is given by / = ig(p — p ), where we 
have set the electron charge e = 1 for convenience. The 

A 

time dependence of the expectation value (I) t can be ob- 
tained by solving Eq. (7) and performing the inverse 
Laplace transformation. For time t — > oc, the result is 

A 2g 2 TuTpB 

[ g 2 r D B + [ g *B + r D + 2 7 r£ + ( 7 2 + n 2 W%] 

b = 7 + ( 7 2 + n 2 )r D , (8) 

where g 2 tt and g 2 j are the exciton frequency shift and 
decay rate, respectively. The derivation of the current 
equation is closely analogous to the spontaneous emis- 
sion of phonons in double dots|2^], in which the correla- 
tion functions /XfXj,^ is given by the elect ron-phonon 
interaction. 

Since the stationary current through the quantum dot 
depends strongly on the decay rate 7, the results of a 
quantum dot inside a planar microcavity is numerically 
displayed in Fig. 2. In plotting the figure, the current 
is in terms of 100 pA, and the cavity length is in units 
of Ao/2, where Ao is the wavelength of the emitted pho- 
ton. Furthermore, the tunneling rates, Tu and Tp, are 
assumed to be equal to 0.2 7o and 7o , respectively. Here, 
a value of 1/1. 3ns for the free-space quantum dot decay 
rate 7o is used in our calculations jl9| . Also, the planar 
microcavity has a Lorentzian broadening at each reso- 
nant modes (with broadening width equals to 1% of each 
resonant mode) jl7| . As the cavity length is less than half 
of the wavelength of the emitted photon, the stationary 
current is inhibited. This is because the energy of the 
photon generated by the quantum dot is less than the 
cut-off frequency of the planar microcavity. Moreover, 
the current is increased whenever the cavity length is 
equal to multiple half wavelength of the emitted photon. 
It represents as the cavity length exceeds some multiple 
wavelength, it opens up another decay channel abruptly 
for the quantum dot exciton, and turns out that the cur- 
rent is increased. With the increasing of cavity length, 
the stationary current becomes less affected by the cavity 
and gradually approaches to free space limit. 

To understand the inhibited current thoroughly, we 
now fix the cavity length equal to Ao/2 and vary the 
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FIG. 2: Stationary tunnel current, Eq. (8), as a function of 
cavity length L c . The vertical and horizontal units are 100 
pA and Ao, respectively. 

exciton energy gap, while the planar microcavity is now 
assumed to be perfect. The vertical and horizontal units 
in Fig. 3 are 100 pA and 2/ic/Ao, respectively. Here, Ao 
is the wavelength of the photon emitted by the quantum 
dot exciton in free space. Once again, we observe the 
suppressed current as the exciton energy gap is tuned 
below the cut-off frequency. The plateau features in Fig. 
3 also comes from the abruptly opened decay channels for 
the quantum dot exciton. From the experimental point 
of view, it is not possible to tune either the cavity length 
or the energy gap for such a wide range. A possible way 
is to vary the exciton gap around the first discontinu- 
ous point 2/ic/Ao. Since the discontinuities should smear 
out for the real microcavity, it is likely to have a peak 

A 

if one measures the differential conductance d(I)/de as a 
function of energy gap e. 
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FIG. 3: Stationary current as a function of exciton energy 
gap e. The cavity length is fixed to Ao/2. The current is in 
units of 100 pA, while the energy gap is terms of 2/ic/Ao. 

Our proposal can also be used to measure the super- 
radiance of quantum dots in an electrical way. Consider 
now the system containing two quantum dots with a dis- 
tance d as shown in Fig. 4(a). One of the obstacles in 
measuring the superradiance between the quantum dots 
comes from the random size of the dots which result in a 
random distribution of energy gap. This can be overcome 
by growing two gates above the quantum dots. The en- 
ergy gap and the orientation of the dipole moments can 
be controlled well. Analogous to the two- ion system [p3[, 



= -i 9 -(p{z) -p*(z)) + 2-^(l/« - n v (z) - n D (z)) 

= l( p (z)-p*(z))-2^n D (z) (7) 

z z 

= -ig{nu(z)C £ (z) - n D (z)C*_ £ {z)} - T D p(z)C e {z) 
= ig{ nu {z)C:{z) - n D {z)C-Jz)\ - T D p* (z)C* (z) . 



FIG. 4: (a) Proposed device structure. Two InAs quantum 
dots are embedded in a p-i-n junction. Above the quantum 
dots are two gates, which control the energy gap and orien- 
tation of the diploes. (b) Stationary current (in units of 100 
pA) through the superradiant (solid channel) and subradiant 
(dashed line) channel as a function of dot distance d (in units 
of A ). 



In conclusion, we have proposed a method of detecting 
the Purcell effect in a semiconductor quantum dot sys- 
tem. By incorporating the InAs quantum dot between a 
p-i-n junction surrounded by a planar microcavity, the 
Purcell effect on stationary tunnel current can be exam- 
ined either by changing the cavity length or by varying 
the exciton energy gap. Second, it is also possible to ob- 
serve the superradiant effects between two dots by using 
present model. The interference features are pointed out 
and may be observable in a suitably designed experiment. 
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